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Abstract 

In this paper we consider excited state g~functions, that is, overlaps between 
boundary states and excited states in boundary conformal field theory. We find 
a new method to calculate these overlaps numerically using a variation of the 
truncated conformal space approach. We apply this method to the Lee- Yang 
model for which the unique boundary perturbation is integrable and for which 
the TEA system describing the boundary overlaps is known. Using the truncated 
conformal space approach we obtain numerical results for the ground state and the 
first three excited states which are in excellent agreement with the TEA results. 
As a special case we can calculate the standard g-function which is the overlap 
with the ground state and find that our new method is considerably more accurate 
than the original method employed by Dorey et al. 
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1 Introduction 



If a conformal field theory is defined on a domain with a boundary then boundary conditions 
need to be defined. In the simplest case of the unit disk, the boundary conditions are im- 
plemented through a boundary state. Such boundary states are not normalisable, but are 
uniquely determined by their overlaps with the normalisable, finite-energy states of the bulk 
system. The exact forms of the boundary states of all the conformal boundary conditions for 
many conformal field theories, and hence all these overlaps, have been known for some long 
time. 

In this paper we consider non-conformal boundary conditions generated by a perturbation 
of a conformal boundary condition by the integral of a local field along the boundary. The 
resulting perturbed boundary condition can also be described by a boundary state which is 
uniquely specified by its overlaps with the bulk states. If the perturbation is integrable then 
it may be possible to calculate the overlaps using the Thermodynamic Bethe Ansatz (TBA). 
For such systems there is a distinguished basis of states and for each state there is an integral 
equation for the overlap. If the perturbation is not integrable then one needs an alternative 
method to calculate the overlaps and determine the boundary state. 

One important numerical method to study perturbed conformal field theory is the Trun- 
cated Conformal Space Approach (TCSA). This was initiated by Yurov and Zamolodchikov 
1^] and used to study the finite-size spectrum of bulk perturbations of conformal field theo- 
ries. It was adapted by Dorey et al to study the finite-size spectrum of boundary-perturbed 
conformal field theory and it was shown by Dorey et al in how this could be used 
to calculate the g-function, or ground state overlap with the boundary state. It was not 
possible, however, to use this method or a modification to calculate overlaps of the boundary 
state with excited bulk states. 

Here, we present an alternative method based on the TCSA inspired by the work of 
Zamolodchikov on the partition function of a perturbed conformal field theory on a sphere 
1^] where the partition function on a sphere was calculated using the Schrodinger equation 
for a truncated system. The sphere was conformally mapped into the complex plane and the 
resulting system was formally the same as the usual truncated conformal space approach with 
a space-dependent coupling. 

We calculate the boundary state overlaps as one point functions of bulk fields on a disk. 
The disk can be conformally mapped to the upper half plane and the one-point function of a 
bulk field is mapped to the expectation value of the bulk field in a state corresponding to a 
half-disk with a time-dependent coupling to a boundary field. We then truncate the Hilbert 
space and numerically integrate the perturbation to find the state for the half-disk. Once the 
half-disk state has been found, it is only necessary to calculate the matrix elements of the 
required bulk field to find the overlap of the bulk state with the boundary state. 

There are some technical difficulties applying this method to higher excited states arising 
from the form of the TCSA expression for the bulk field expectation value which is similar 
to a Fourier expansion in the argument of the insertion point. Excited states correspond to 
repeated derivatives of the field and the derivatives of the Fourier series do not themselves 
converge, this can be handled by forcing convergence of the Fourier series (at the loss of some 
information) which gives numerically satisfactory results at the expense of some complication 
in the calculational scheme. 

We apply this method to the simplest perturbed boundary conformal field theory, the Lee- 
Yang model, which was the model studied in ^ and has the advantage that the overlaps 
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can also be calculated using the TBA method and we find excellent agreement. 

The layout of the paper is as follows. In section 2 we outline the general truncated space 
method and in section 3 we apply it to the ground state and first excited state in the Lee- 
Yang model. In section 4 we consider the next two excited states and explain how to force 
convergence of method. We outline the TBA method used in an appendix. 

2 gf-functions and perturbed conformal field theory 

Our aim is to calculate the overlap of a bulk state with the boundary state corresponding to 
a perturbation of a conformal boundary condition on the boundary of the disc \z\ = Rhy 
the insertion of exp(— A J ■0(x)dx) on this boundary. We take V to be a primary boundary 
field of conformal dimension h. We shall denote such an overlap by Qi where i labels the 
particular bulk state. The most natural way to calculate this one point function is to expand 
the exponential and consider it as a perturbative expansion with repeated insertions of the 
boundary field. If /i < 1/2, then such a function Qi will have a regular expansion in powers 
of {XRy) for small A (where y = 1 — his the mass dimension of the coupling constant A), but 
for large A the dominant term is non-perturbative: 

\ogGi ~ -RfB + ... , (2.1) 

where fs is the boundary free energy per unit length and R = 2ttR is the length of the 
boundary. While this is the natural consequence of a perturbative expansion, it is many ways 
more natural to consider the state perturbed boundary to interpolate between two conformal 
boundary conditions which requires subtracting the free energy term; this gives the definition 
of the (non-perturbative) g-function gi as 

g,{\Ry) = exp(2^/iji?) GiiXRy) . (2.2) 

We shall define the perturbed boundary state {B{\)\ by the conformal normalisation so that 
the overlap gives with the bulk state \(pi) the perturbative ^-function, 

gi{XRy) = {B{\)\^i) . (2.3) 

The overlap {B{X)\(pi) can be realised as the insertion of the field (/?i(0) in the disk with 
boundary condition B{X), where \ip)^ = ^i{0)\ 0) and it is this definition which is the basis of 
the TCSA method described below. 

2.1 The TCSA approach to the disk 

In order to use the TCSA method, we need to find a way in which the boundary field can 
be constructed as an operator acting on a Hilbert space (which can be truncated), or equiva- 
lently in which the system has a Hamiltonian description where the boundary field enters the 
Hamiltonian. We first outline the operator realisation and then the Hamiltonian description. 

To have an operator action, once must choose "equal-time" surfaces to which one can 
associate a Hilbert space of states on these surfaces. For the boundary field to act, these 
surfaces must have boundaries, and the most natural surface is a line segment corresponding 
to equal time slices across an infinite strip. This is conformally equivalent to the upper half 
plane with the equal-time surfaces being semi-circles of constant radius, known as radial 
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Figure 1: The maps i : z ^ w and t : w ^ showing the "equal-time" quantisation surfaces 
as dashed hues. 



quantisation. We can map the interior of the unit disk (with coordinate z) to the upper half 
plane (with coordinate w) using the map 



i : z ^ w 



z+l 



(2.4) 



This map is shown in figure || The one-point function in the perturbed boundary condition 
(v(0))a becomes the matrix element of a field between states representing the interior of 
the unit semi-circle (with perturbed boundary condition) and the exterior (with perturbed 
boundary condition). Using the transformation property of primary spinless boundary fields, 



L : miz] I— > 



2\ h 



{w) , (I){exp{i9)) ^ (^^— J Ht) , t = tan(0/2) , 
we can write the state corresponding to the interior of the semi-circle of radius r as 
\^{r)) = Pexp ( - A j\<t>{t) + ct>{-t)){j^)y dt) | 0) , 



where y = 1 — h and P denotes radial ordering. The overlap {B{X)\(p) is then 

(S(A)|^)= const. (M/(1)|<(^(0))|^(1)), 



(2.5) 



(2.6) 



(2.7) 



where (const.) is a constant normalisation factor and where L{(p{0)) is the image of the field 
(p{0) under the map t; if f{z, z) is a spinless quasi-primary field of total conformal weight A, 
this is the simple rescaling 



Lifiz,z)) 



dw 



dz 



(2.^ 



otherwise it is more complicated. 

The state |^(r)) satisfies the differential equation 



dr 



\^(r)) = -X{cPir) + ^i-r))ij^)y\^ir)). 



(2.9) 



and we can use ( p.9| ) as the basis for a TCSA calculation. We can truncate the Hilbert space 
of the system at excitation level N and solve the first order matrix differential equation for 
|^(r))jy using the finite matrix truncations of {(j){r) + (j){—r)). We can then take the matrix 
elements of /-(y?) in the truncated state |^'(l))jy and this is the TCSA approximation to the 
overlap. 
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2.2 Realisation as a Hamiltonian system 



The calculation of the previous section can be recast in terms of a perturbed Hamiltonian (in 
the spirit of Q) as follows. It is easiest to see the relation by mapping the disk to an infinite 
strip, by the further map 

T : w ^ u = log{w) , (2-10) 

which maps the original disk to the strip < Im(u) < vr (see figure |l]. If u = x -\- iy then the 
lines X = const, can be taken as equal time slices and the Hamiltonian propagates in x along 
the strip. The equation ( ^.9| ) is transformed into the Schrodinger equation in the interaction 
picture, 

^\^{x)) = -Hj{x)\^{x)) , (2.11) 

where Hj{x) is the x-dependent Hamiltonian 

Hi{x) = Asech(x)^((/)(x) + 0(x + in)) . (2.12) 

This problem can be also be formulated in the Schrodinger picture in which case the relevant 
Hamiltonian is 

H{x) = Lo + A sech{x)y{^{0) + (j){i7r)) . (2.13) 

It is this formulation which is directly comparable to the Hamiltonian formulation of the bulk 
model on the sphere in Q. In practice we have used (|2.9| ) to perform our TCSA calculations 
as it is numerically more stable. 

We now turn to the calculation of the ground state and first excited state overlaps in the 
Lee- Yang model. 



3 The Lee- Yang model 

The Lee- Yang model has been used extensively as a test of exact and numerical methods 
because it is in many ways the simplest, and the convergence of the TCSA calculation of its 
spectrum is amongst fastest of all known models. In this context it is the simplest because 
it has only two conformal boundary conditions and a single integrable perturbation which 
connects them. It also has the fewest states at any particular excitation level so that the size 
of any TCSA truncation will be smallest in the Lee- Yang model. The details of the conformal 
and integrable descriptions of the boundary conditions can be found in [||, ^ and we only 
recall the most important details here. 

We start with the description of the theory in the bulk. The Lee- Yang model is the 
minimal conformal field theory with c = —22/5. It is non-unitary and there are two highest 
weight representations of the Virasoro algebra which are relevant, with conformal weights 
— 1/5 and 0. If we denote these two highest weight representations by H_i/^ and Hq then the 
Hilbert space of the bulk theory is 

nbuik = (Ho ® Ho) e (//_i/5 «> ^-1/5) • (3.1) 

We take the bulk highest weight states to be | 0) and \(p) of conformal weight and —2/5 
respectively. The SL{2) invariant state is | 0) but the ground state in the bulk theory is \ip). 
It is convenient to normalise these states as 

(0|0) = -1, {^\^) = l. (3.2) 
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The state \(p) is associated to the primary field (p{z,z) and satisfies 

|^) = ^(0,0)|0) . (3.3) 

The Lee- Yang model has two conformal boundary conditions. We will adopt the conventions 
of in which these are called 1 and $ and the states associated to them on the unit disk 
are 

{B^\ = g^{{^\ - Z$((0| , = g^M - Zi((0| , (3.4) 

where ga is the g-value of the boundary condition a, is the value of the disk partition 
function with boundary condition a and ((^| are the usual Ishibashi states [||. 

As described in [||, there is a single integrable flow between the boundary conditions $ 
and 1, the perturbed boundary conditions being denoted by $(/i) where h is the coupling to 
the unique relevant boundary field on the condition <I>. This field has conformal weight —1/5 
and is denoted by (f). 

The particular boundary conditions we are considering, both perturbed and unperturbed, 
are rotationally invariant and so the only states which have non-zero overlap must have zero 
spin, that is they must have equal Lq and Lq eigenvalues. The lowest lying zero-spin states, 
their conformal weights and the fields to which they correspond are in table I. We'll denote 
the overlap with the i-th excited state as Qi and the associated (/-function as gi. It is a 
property of the conformal boundary states that the overlap with a normalised bulk state only 
depends on the representation in which it lies, so that in particular we have 

ao(0) = ^2(0) = ^3(0) = ^5(0) = {B^\ip)=g^ , gi(0) = ^4(0) = (S$|0) = . (3.5) 



state weight field 

W) -2/5 
|0) 

-|L_iL_i|(^) 8/5 

-f|L_iL_iL_iL_i|(^) 18/5 
-^L_2L_2|0) 4 

-i§L_iL_iL_iL_iZ_iZ_i|(^) 28/5 

Table 1: Lowest lying spin zero states in the Lee- Yang model 

We also need the description of the $ boundary condition in the upper-half plane op- 
erator picture. In this formulation, the unperturbed Hilbert space is that of a strip with 
the <I> boundary condition on both sides and is the direct sum of the two Virasoro algebra 
representations in the Lee- Yang model, 

T~(-strip = Hq(B H^i/^j . (3.6) 

We denote the highest weight states in the UHP picture as |0) and respectively. The 
ground state is \<j)) which corresponds to the insertion of the boundary field (j){0) at the origin 
of the model on the upper half plane. We again choose to normalise these states as 

(0|0) = -1, = 1. (3.7) 



ip{z) 
1 

— I dd(p{z, z) 
-§d^d^^iz,z) 

-'inz)f(-z) 
_i|93av(^,^-) 
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The expectation value of any operator O on the disk with the conformal <I> boundary condition 
is now given by 

{O)^ = -Z^{0\l{O)\0) , (3.8) 

where i is the map (^^) from the disk to the UHP. In particular we recover the partition 
function Z$ for the identity operator = 1, 

(1)^ = -Z$(0|0) = Z$. (3.9) 

For the ground state entropy we use the operator O = (^(0,0) and i{ip{0,0) = 2~'^/^ip{i, -i). 
With the result {0\(p{w,w)\0) = -(*)5^(2r sin 61)2/5 we get 

( (^(0) )^ = -Z$( I 2-2/5^(,, -i)\0) = Z^ = . (3.10) 

We now turn to the overlaps with the excited boundary states. 

The first overlap studied was the ground state overlap, the original (7- function of p| for 
which a TBA formulation was first described in Q. Subsequently, it was shown in Q^that 
the 5f-function was the first of an infinite set of excited state g-functions which arose as the 
eigenvalues of generalised transfer matrices satisfying a so-called T-system. This T-system 
has been put into an orthodox conformal field theory context by Runkel @ . In the case of the 
Lee- Yang model the T-system is identical to the so-called y-system and Bazhanov et al also 
showed how to find the TBA system for excited state y-functions and solved it for the first 
excited state. The ground state y-function was compared to an alternative TCSA calculation 
of Qo in ^ with reasonable agreement. Here we state the expected relation between the Q 
and Y functions: 

Yo{e)=y^{B^{h)\ip) , Y,{9) = -^{B^{h)\0) , (3.11) 
5* Z^ 

where yo = ^{1 + \/5), yi = ^(1 — v^) and we have suppressed the /i-dependence of the state 
1^(1)). The y-functions depend on a rapidity 6 which is related to /i by 

K2.r^ = -1 K e'"' , K = -vr^/^ (r(f )r(|))V2 [y^) = "O-^^^^g... (3.12) 

The extra factor of (27r)^/5 between the formula here and that in Q arises because the 
boundary of the unit disk has length 27r and it is the combination hR^^^ which is relevant to 
this equation. 

We now describe the TCSA approximation and in particular the lowest level truncation 
where the UHP Hilbert space is restricted to just two states. 



3.1 TCSA method and results 



We have used the TCSA approximation to the differential equation ( 2.11 ) to calculate ap- 
proximations to the state |^'(1)) and then calculated the exact expectation value of the fields 
{p{i) and IL in these approximate states for truncation levels up to 20, that is for systems of 
size 2 up to 323. To illustrate the method, we consider here the case of truncation to level 0, 
that is we truncate the Hilbert space to just the two highest weight states | 0) and and 
construct the TCSA equations for the upper half plane state (|2.9|) explicitly. 
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We take the state \^{r)) to be truncated to just the highest weight states \<j)) and | 0), so 
that it is determined by just two functions /^(r) and fo{r), 

\^{r))=U{r)\<p) + /o(r)|0) . (3.13) 

The actions of the field (/)(r) on these states are (for r both positive and negative) 

4>{r)\0) = |(/<)+rL_i|0) + ... , (3.14) 
0(r)|c/.) = -|r|2/5(|0) + ^r2L_2|0) + ...) 

+ c<^/|r|i/5(|</,) + lrL_i|,/.) + ...) . (3.15) 

This means that truncating (^]^) to the two states gives the differential equation 

d;^ U(0 j " " \2 2c^/rV5j \f^ir)) ' ^^'^^^ 
which should be solved subject to the initial conditions 

/o(0) = 1 , /^(O) = . (3.17) 

The result is the TCSA estimate 

1^(1)) = a|0) + 6|</.) , a = /o(l), 6 = /i(l). (3.18) 

From this we can calculate the TCSA estimates of the two ^-functions, 

Qo = 2-2/5 (^(1)1 (^(i,-i) 1^(1)) , (3.19) 

Gi = (^(1)1^(1)) . (3.20) 

To complete the calculation of the Y-functions, we need the one-point functions of the field 
ip{w,iv) on the UHP, which are in appendix^ and the re-scalings from ( |3.11| ), resulting in 



'•"H^^, I g-9/5 (^^-i)r(i^)r(i) 



= 1.618034..(a^ - 2.50895.. a5 + 1.50258.. 6^) , (3.21) 
Yi = yi (a^ - 6^) = -0.618034.. (a^ - 6^) . (3.22) 

We used Mathematica [^] built-in routines to solve the matrix differential equations with 
satisfactory results. Already the accuracy obtained using just a truncation to 2 states is 
quite surprising. We list some numerical results in table |2| and in figure ^ we plot the 2-state 
truncation estimates for the functions Yq and Yi together with the exact results calculated 
using the TBA equations of |p and the results of second and third order perturbation theory 
(see appendix^. 

Although the agreement with the functions Yq and Yi as plotted in figure 1. of |P is 
impressive, this is in main due to the dominance of the functions Yi by the boundary free 
energy term which means that the Y functions grow rapidly and are hard to show on this 
particular plot. To remove this factor, we instead plot the the (7-functions gQ{hRy) and 
gi{hRy) which interpolate the UV and IR values of Qa and Za respectively, using 

2TTfBR = -TT\2h/hc\^/^ . (3.23) 



8 



log/i 


a 


b 


VTCSA 


vexact 


yTCSA 


■y^ exact 


— oo 


1 





1.61803 


-0.618034 


-6 


0.999783 


-0.0237672 


1.65316 


1.65471 


-0.617952 


-0.617953 


-5 


0.999783 


-0.0237672 


1.71517 


1.71948 


-0.617417 


-0.617433 


-4 


0.998362 


-0.0669088 


1.89480 


1.90719 


-0.613244 


-0.613387 


-3 


0.987093 


-0.200034 


2.47539 


2.51448 


-0.577453 


-0.578895 


-2 


0.886238 


-0.704015 


5.00870 


5.16228 


-0.179094 


-0.201836 


-1 


-0.379141 


-3.85251 


30.3869 


31.0587 


9.08392 


8.03592 



Table 2: Numerical results for the two-state TCSA approximation 



N 


90 


91 


92 


93 


94 


95 





0.853057 


0.0387796 


0.430042 


0.561775 


0.106314 


0.687637 


2 


0.889914 


0.0422042 


0.455962 


0.653033 


0.103007 


0.993606 


4 


0.878654 


0.0428922 


0.437669 


0.487325 


0.125051 


0.688144 


6 


0.882040 


0.0431683 


0.451097 


0.745018 


0.076786 


2.72391 


8 


0.879529 


0.0433135 


0.440139 


0.400679 


0.153111 


- 0.385525 


10 


0.880638 


0.0434018 


0.449176 


0.830633 


0.048286 


6.80768 


12 


0.879593 


0.0434606 


0.441394 


0.322728 


0.181909 


- 3.52987 


14 


0.880125 


0.0435024 


0.448145 


0.907610 


0.019363 


14.0118 


16 


0.879565 


0.0435335 


0.442140 


0.251577 


0.210946 


- 9.55319 


18 


0.879871 


0.0435574 


0.447506 


0.978097 


-0.009753 


25.0185 


20 


0.879526 


0.0435764 


0.442631 


0.185561 


0.240248 


-19.1595 




0.8794 


0.04368 


0.4449 


0.6020 


0.1079 


- 6.93 


ooiQ) 


0.8794 


0.04358 


0.4449 


0.6143 


0.1008 


0.6912 


TBA 


0.879214 


0.0437266 


0.444889 


0.614126 


0.100843 


0.690846 



Table 3: The values of go, 91,92,93, 94 and (75 for log h = —2 estimated from TCSA at trunca- 
tion levels < < 20, the values obtained by Pade extrapolation of the series (P) and (Q), 
as well as the exact value from TBA. 



In figure ^ we show tha TCSA estimates (for truncation levels and 20) and TBA values of 
qq and gi, along with the result of a Pade approximant extrapolation in truncation level to 
N = 00. As can be seen in figure 3(a) , the estimate of go and gi from the level zero truncation 
is already quite good and convergence improves rapidly with increasing truncation level, so 
that for level 20 it is qualitatively accurate up to log/i ~ —0.5. The main error seems to come 
from errors in the extensive boundary energy per unit length which can be attributed to an 
effective coupling constant renormalisation in the TCSA method [^]. Rather than attempt 
to incorporate these effects in the differential equation, we can remove the error from the 
boundary free energy by looking at the ratio of the (7-functions, which is plotted in figure 

We can attempt to improve the accuracy of our estimates by using an extrapolation 
method to estimate the N ^ 00 limit of the TCSA data. We have used two ways of approach- 
ing the — 7- 00 (described in section ^) limit and three different extrapolation methods. In 
all cases, all three extrapolation methods gave the similar results, with Wynn's method |11| 
based on Pade approximants being the best. The two different ways of approaching the limit, 
while very similar for gQ and gi, differed considerably for some higher excited states. The 
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Figure 2: The two lowest lying y- functions in the Lee- Yang model, as calculated from TBA 
equations (solid lines), and the TCSA approximations calculated using the space truncated 
to 2 states (points) plotted against log(/i(27r)^/^) for comparison with [^. Also shown are the 
second order perturbation theory expressions (dashed lines), third order (dot-dashed) and 
sixth order (dotted) - see appendix ^ for details. 



simplest method, considering the level approximation to be given by the full TCSA ex- 
pression (^^), failed to be effective for fields of weight greater than 2, and we had to use a 
second approximation to get an effective extrapolation to = oo. We illustrate this in table 
^ with the TCSA estimates of qq and gi at log(/i) = —2 for < < 20 together with the two 
Pade extrapolations labelled oo^^^ and 00^*^^ and the TBA value. We also include the values 
of 92, 93, OA and 55 which will be discussed in the next section, together with a discussion of 
the two extrapolations. 

To conclude this section, we discuss the comparison with the previous TCSA estimation of 
(7o in Q. In that paper, was estimated from the scaling with the strip width of the cylinder 
partition function calculated as a trace over the space of states on the strip. This was a rather 
convoluted calculation which nevertheless showed good agreement with the exact functions. 
In figure 3(a) we include the old TCSA estimate from the level 18 truncation (using the exact 
boundary free energy per unit length as shown in [^); as can be seen, the new method is much 
better. The old TCSA estimate could be improved by fitting the boundary free energy term 
to account for the renormalisation of the coupling constant, as can the new TCSA estimate 
by extrapolation in N , and the new method again gives better results. 



4 Higher excited states in the Lee- Yang model 

We can continue to apply the same method to the higher excited states in the Lee- Yang 
model. The next four states are given in table |l] and three of the corresponding fields are 
given by simple derivatives of the primary field {p{z,z). We can calculate these by directly 
taking the derivative of the one-point function on the perturbed disc of (p{z,z), which we 
discuss in appendix |^. The fourth excited state corresponds to the field T{z)T{z) and we 
discuss its matrix elements in appendix The results for the excited state functions 52 to 
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(a) The function go (b) The functions gi 

Figure 3: The TCSA estimates for the g(- functions go and gi level (open squares), 20 (open 
circles), and extrapolated (solid circles), compared with the exact results from TBA (solid 
lines). Also shown are the UV and IR values as straight lines and the level 18 result for gQ 
from the method of (dashed line). 

0.5 

0.0 
-0.5 
-1.0 



-5 5 

Figure 4: The TCSA estimate of the ratio gi/go from truncation at level 20, together with 
the exact results from TBA (solid lines). Also shown are the UV and IR values as straight 
lines. 

(75 at log/i = —2 in table ^. As can be seen, these are convergent for the first three states 
but oscillatory and divergent for g^ to g^. The behaviour of these divergent cases is of an 
increasing oscillating series, typical of a divergent perturbation expansion. 

This immediately suggests that one might be able to re-sum the divergent series and 
obtain useful results. We have tried Euler, Cesaro and Fade extrapolation (using Wynn's 
method [^]) and these all gave similar results. They all improve the results ofr go to 52 > 
are all able to produce a result correct to a few percent for the divergent series 53 and 54 
but all fail on g^. As a consequence we considered different ways to approach — )• 00 limit, 
that is we constructed a different TCSA approximation at each level which would agree for 
convergent series, but could give better results for divergent series. Looking at the matrix 
elements {a\Oi\f3) of the fields Oi it seemed that these were organised by the total excitation 
level of the states \a) and |/3) rather than the maximum level. We consequently defined two 
series as follows. We first denote the contributions at excitation level n to the level N TCSA 
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approximation \^n) to the state |^(1)) by \^n), so that 

AT 

\^N)=^\i^n), (4.1) 
n=0 

For any operator O, we then define two series of polynomials 

TV 

m=0,n=0 
N 

QnW = E A"^+"(^^|0|V™> . (4.2) 

m+n=0 

The direct TCSA approximation at level is given by P^{1). If the TCSA approximation 
converges m N then both P^(l) and Q%{1) converge to the same answer; if the TCSA 
expression diverges then the two polynomials can differ considerably and we can take the 
Fade approximants to Pat (A) and QnW at A = 1 and see if they provide useful estimates of 
the answer. As an example, we consider O = (— 11/5)T(0)T(0) and the truncation to level 
12 at h = 1. The polynomials are 

P[^T = 657.76 - 1077.2 A^ + 3234.6 A'^ - 5513.1 A^ + 8256.0 A^ - 11183 A^^ + 14552 A^^ , 
Ql^ = 657.76 - 3671.3 A^ + 8657.4 A^ - 14697 A*' + 21200 A^ - 27983 A^° + 35138 A^^ , 

Pj^^(l) = 8926.8 , Qi2il) = 19301. . (4.3) 

We see that neither of these is close to the exact result of -113.6. Instead we take the Fade 
approximants 

657.76 - 234.64 A^ + 1388.4 A'' - 1277.3 A*^ 

1 + 1.2810 A2 - 0.70901 A^ - 1.0206 A^ ' 
657.76 -2011.4 A2 + 766.95 A^ - 127.37 A^ 
1 + 2.5237 A2 + 2.0903 A^ + 0.60087 A^ ' 

and find the Fade approximant estimates and the exact (TBA) answer are 

pff (1) = 968.9 , 9^2^(1) = -114.9 , GJ^^ = -113.6 . (4.4) 

We see that the Q-approximation is reliable and the P-approximation is not. This is repeated 
for the other divergent cases of and 5(4 - the Q-estimate is always better than the P-estimate. 
It is worth noting that for the convergent cases there is very little difference between the two, 
but that the P-estimate is very slightly better than the Q-estimate since the P-polynomial 
includes more (convergent) matrix elements than the Q-polynomial. 

We show the resulting plots for the ratios Qi/go of the Q-extrapolated TCSA results at 
level 20 in figure |. We see that the ratios, which are the physical content of the 5- functions, 
are all well described by the Q-extrapolated TCSA data up to approximately the same value 
of log h ~ 2.5, which is the value at which the TCSA Hamiltonian on a strip also fails to give 
an accurate description for the 5th excited state of the spectrum. 



ql2^W = 



12 



-5 







5 



Figure 5: The extrapolations of the TCSA estimates for gi/go plotted against log(/i) for 
i = 1,2,3,4 and 5 (labelled by i). Also shown are the TBA values (solid lines) and the only 
allowed values in the IR and UV (straight lines). 



5 Conclusions 

We have found a new and effective way to calculate the overlaps of perturbed boundary 
states with bulk states giving excited state functions, or equivalently T- functions. We have 
studied these in some detail in the Lee- Yang model (in which the T-functions are identically 
equal to the y-functions) and shown that the estimates are good for the whole range of 
values for which the TCSA method on the strip is applicable. Even using a truncation to just 
two states, the results are better than second order perturbation theory. It is disappointing, 
although perhaps not surprising, that the TCSA estimates computed directly are divergent 
for the states of conformal weight greater than 2, but we have found an efficient extrapolation 
method which works effectively for all the states we have considered. 

This new TCSA method can be applied to any boundary perturbation for which the 
boundary couplings are known and is not restricted to integrable models, as was the case 
here. The method also has immediate scope for generalisation to include massive, or bulk, 
perturbations, which should help check any conjectured TBA equations for the bulk excited 
state (^-functions which are complicated by the extra infinite series required to account for the 
bulk corrections 0]. It can also be applied to study defect matrix elements (which are 
boundary state overlaps in the double model through the usual mirror trick) and calculate 
the reflection and transmission coefficients defined in [|l^ ] and help clarify the nature of the 
endpoints of the defect flows studied in []l9| . 
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A The Lee- Yang conformal data 



The coefficients appearing in the boundary states (3^) are 



5$ 



V5 
V5 + 2 



V5 



1/4 



1/4 



0.570017... , Z 



1 



\/5 + l 



2^/5 



1/4 



1.17319... , <7i 



2V5 



1/4 



= -0.922307... , (A.5) 
0.725073... . (A.6) 



The various structure constants used in this paper are 



5 + V5 



1/2 



• a 



2.96585.. , 



l + \/5 



2 

l + \/5 



1/2 

• a = -1.98338.. 



3/2 



(A.7) 



2.05817.. 



where 



a 



r(i/5)r(6/5) 



1/2 



1.55924.. . 



r(3/5)r(4/5) ~ (A-^ 

Finally, the one-point functions of ip{w,w) on the UHP are given in terms of r = \w\ and 
9 = aig{w) as 



-_ -W^i (2r sin 0)2/5 

-_ -(*)s^r2/5(2r sin 0)1/5 
= (2rsin0)V5 

= -cr2/5(2sin0)i/5F(cos2(0)) 

^ viO' 10' 2^-^)- 



(A.9) 



{0\ip{w,w)\0) 

{ I w)\(j)) 

{ (j) I ^) I ) 
{(l)\ip{w,w)\(l)) 

where c = ^^^^^^^ and F(x) 

B Perturbation theory 

It is straightforward to obtain one or two orders of the perturbation epansions of the functions 
Qi. For the 5- function Qq this was given in |p and the conformal theory perturbation discussed 
in [Q. Here we write them in terms of hR^/^ rather than h{2'KR)^/^: 

Go = 5$ (1 + aihR'^/^ + a2{hR'^/^f + 

1/2 



ai = (27r) 



^/5r(2/5)r(6/5) 



012 = OL 



2/^/5. 



2cos(^/5)r(4/5)2 

The function Qi, the disk partition function, can be expanded out to third order as 



Gi = zJl - UhRyf {2tt) 



22^^^^(^/l«) +i(.i?^)3(2.)3^(^^/^0) 



V^r(6/5) 



r(ii/io)3 



+ 



(B.IO) 



(B.ll) 
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Similarly, G2,G3 and G5 can be expanded to second order and Ga to third order based on 
conformal perturbation theory; beyond these orders the integrals become difficult to evaluate 
exactly, although use of the Y-function relation ( [E.30D does allow one to calculate one further 
order for Gq and three further orders for Gi as the coefficients in the perturbation expansion 
are severely constrained by this relation and so obtain 

Go = 5*(l + «ia; + ^(ai)V + ^^|^(ai)V + ... ) (B.12) 



Gl = Z$(l + /32x"+/33X^-^X^ + ^^/32/53X 

(l + V5)(/32)^ + 2^/5(/33)^ 
10 



+ (B.13) 



where x = hR^^^. These are also shown in figure ^ 

It is rather harder to calculate the z-dependent function ( ip{z,z) ) as the integrals again 
become harder. To first order we have (for a disk of radius 1) 

{^{z,z))=Z^{l-\z\r^' + (2vr)/.(*)i?^(l-|.|2)-VSFi(-l,-l;l;|zp) + .. 

(B.14) 



C ip{z, z) expectation values 

Of the first four excited states given in table (I), thr ee are derivatives of the primary field 
(p{z,z). We can calculate the overlaps with these states by directly taking the derivative of 
the one-point function of {p{z,z). Here we show how this is defined in our TCSA approach. 
We use the coordinates z on the unit disk and w = rexp(i^) on the UHP, related by ( p. 41 ). 
In these coordinates 



1 — + 2ir cos ( 
1 + r"^ + 2rsm6 



1 + — 2r sin 6 
1 + r2 + 2r sin 9 



dw 



dz 



i(l + r2 + 2rsin6l) . (C.15) 



The general form of the expectation value of an operator O on the disk is given by and 
the transformation property of the primary field ip{z, z) is given by (|2.5D, so that 



( ^{z,z) 



/disk 



;i(l + r2 + 2rsine)) ( ^{w,w) 

-2/5 



'UHP 



Z<s> (i(l + r2 + 2rsin^)) (^(l)|(^(re^^ re-^^)|^(l)) . (C.16) 



We then used 



to arrive at 



2/5„Lo+Lo 



ip{e 



ie 



(C.17) 



( ip{z,z) )^,^^ = -Z^ {l{r + l/r + 2sm0)) (vI/(l)|r^°+^o</.(e*^ e"^^^) r-^o-^ol^-a)) . 

(C.18) 

The matrix elements of ip{e^^,e can be calculated in the usual way from the highest weight 
matrix elements (|A.9D . Given these, the TCSA estimate takes the form 



{^{z,z))i% = {l-\z\y' k^ + af(sin 



2\2/5 



-1/5 



N . 



(C.19) 
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where 



d 



/2W =sin(0)-V5 F(i/io, 3/10; 1/2; cos^^) ^ j^^^^ = ^^.^^gyi/^ _^^^-^^g^2/5 (C.20) 



and the coefficient functions have the form 

2N 



2N 



«f = 5;af'^ (r)cos(2j0) , /3f = J] /3f'^ (r) sin(2j0) 
3=0 j=o 

For the level 0, two-state truncation with |'I'(1)) = a| 0) + b\(j)), these functions are 



a. 



The map z ^ w is conformal and the Laplacian in these coordinates is 



dd 



+ 2rsin0 + 1) / 9^ 19 I d'^ 



16 



Q^2 J. Qj. 



89^ 



(C.21) 



(C.22) 



(C.23) 



While the functions converge in A^, the repeated action of dd leads to a failure of conver- 
gence, which in this case is happens first for with two applications of dd. 



D TT expectation values 

The field TT is quasi-primary and transforms under the Mobius map as 

T{z)f{z) ^ \w'{z)\^ T{w)f{w) , 
which gives the expectation value on the perturbed disk as 

' I Jr r"^ + 2r sm{ey ^ 



(D.24) 



( T{z)T{z) ) 



disk 



2 

1 + +2r sin(6')\ 



( Tire^Tire-^'^) 
4 

(^'(l)|T(re^^)r(re-^^)|'f(l)) (D.25) 



The product T{re^^)T{re is regular and its matrix elements can be worked out in the 
truncated conformal space using the normal ordering relation identity appropriate for the 
TCSA method, 



T{z)T{w) = T<o{z)T{w)+T{w)T>o{z) 



+ 



c/2 



{z-wY 



+ 



{z — w)^ z"^ 



T{w) 



w 



z'^{z — w) 



T'{w) , (D.26) 



where T<q{z) = Y.m<o^rnZ ™ ^ and T>o(z) = Y,m>o^rnZ ™ ^. This choice of ordering 
means that the matrix elements of T(e*^)r(e~*^) in the truncated space can be calculated 



exactly using the expansion of the right-hand-side of ( D.26 ) in the truncated conformal space. 
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To work out the level zero TCSA approximation, we just need the expectation value in a 
highest weight state, 



{h\T{z)T{w)\h) 
{h\T{re'^)T{r~'^)\h) 



c/2 



+ 



2h 



+ 



{z — w)'^ zw{z — wY z^nfl 



(D.27) 



c/2 



(2r sin6')2 



m(l-2fe) ^ ^ Cos^g ^o»,x-..; ^ COs"(0)f 



^2 ^16fe(l-4fe) 



Putting this all together, we find the level zero TCSA approximation to the disk one-point 
function of the normalised field is 



^( \nz)fiz) )S:r=(i-N 



+ 6^ — (1 - 72 cos^ + 16 cos^ Q) 
55 



(D.28) 



As is immediately obvious, this does not have the expected rotationally-invariant form 
since the TCSA method breaks rotational invariance. 



E TBA equations in the Lee- Yang model 

The scaling Lee- Yang model is the perturbation of the conformal minimal model A^2,5 by the 
relevant field 

-4sLYM = -42,5 + y <^zAz^\,l[z^ z) 

It has a single particle of mass m in the spectrum with the two-particle S matrix 

sinh 6 -\- i sin % 
S- 3 



sinhf 



I sm 



2n 



The ground state TBA of the SLYM has the form 0] 

/oo jn/ 
— (/^(0-0')log(l + e-^(^')) 
-oo 27r 



TiTBA 
^0 



where 



/oo ^0 
— m cosh 6 log(l + e~ 
-oo 2vr 

\/3 cosh 6 



coslr 6* - 1/4 



The excited state levels can be determined by solving the equations [|13 

k 



e{9) 



ml cosh(0) + log g|g _ j 



/oo J/l/ 
^^(0-e')log(l + e-^(^')) 
-oo 2vr 

1 

i m (sinh 0j — sinh 0*) — / — m cosh 6 log(l + e~ 

.•1 J -oo ^TT 



(E.29) 
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where k is the number of particles present in the state with the singularity positions 9j 
corresponding to their rapidities. For large volumes, the integral terms can be dropped and 
a straightforward calculation shows that the singularity positions have the form 

where the particle rapidities 9j solve the Bethe-Yang equations 

mL sinh 6* j + ^ 5 {Oj - 4) = 27rl. 



where S{9) = -a'^^'^^ 

with Ij e Z (for k odd) or e Z + i (for k even) labelling the particular state in question. 
The Y function of the given level can be defined as 

Y{9) = e^(^) 

and it satisfies the functional relations [^] 

Y(9 + i^)Y(^9-z^)=l + Y (9) (E.30) 

These equations can be solved by iteration: the Bethe-Yang equations give a starting 
point for the position of the singularities, which can then be used to iterate the system 
( E.2E ). For a given position of singularities the function e can be updated by iterating the 
integral equations, and then the positions of the singularities can be refined by solving the 
equations 

3f?e (e-^(^^ ) + 1) = 9m (e'^^^^) + 1) = 

for the 2k real variables given by the real and imaginary parts of 9j. This can be accomplished 
by the multidimensional Newton method; the 2k x 2k derivative matrix of the equations with 
respect to the variables can be expressed by taking the derivative of the integral equation in 

(ill)- 

To describe the states in the minimal model, we consider the so-called kink limit of the 
above TBA system describes a chiral (right moving) version of the c = —22/5 Lee- Yang 
minimal model. The variable 9 must be redefined by the shift 

9^9-log^ 
mL 

Taking the mL — J- limit one arrives at the TBA system 

e{9) = e' + y log ^l^-^j] _ r ^^(0 _ 0') iog(l + e'^^')) 
^ ^ fri 3(9-9*) y_oo 27r^^ ^ ^ 

E{L) = (E.31) 

where the only remaining 9j are those for which the limit 

2 

9j = lim 9j — log 

L-s>oo mL 
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is finite, and the effective central cliarge of the state can be expressed as 



ce. = ^ E (e^"^ - + ^ II dOe^ log(l + e-(^)) 



24A 



where A is the dimension of the operator creating the state and the central charge is c 
-22/5. 

The massless ground state y-function 

Yoi9) = eoie) 

can be found simply by iterating 

/oo jn/ 
—^{e-d')log{l + e- 



which converges very fast. The ground state is created by the operator ip with weight —1/5, 
with Ceff = 2/5. 

For a reliable numerical computation of the massless excited state y-function we adopted 
the following method. We solved the massive TBA equation for a left-right symmetric con- 
figurationQ, and decreased mL to a small value (down to 10^'^) to get a starting position for 
the pseudo energy e and the singularity positions 6j in the massless limit. Then the solution 
was refined by iterating the massless TBA system ( [E.31 ), and the consistency of the solution 
verified by matching the value of Ccff against the value expected from CFT. For the case 
of the second, third and fifth excited states it is only necessary to consider a single pair of 
singularities 9, 9* whose position is 



State 


Ceff 






118 

5 


1.84285939 + 0.006848815 i 


lUv 


238 
3^8 


2.53410313 + 3.03963559 x 10~^i 




5 


2.93773630 + 1.3531641 x lO^^^i 



(up to the accuracy shown). 

The case when the infrared (large mL) configuration contains a stationary particle is 
somewhat exceptional; this is the case for the first excited state which corresponds to the 
identity operator. For large mL the position of the singularity is given by 



1 + ^e-V3mL/2\ 

6 / 



but for small mL the imaginary part of the position of the singularities approaches the values 
i7r/3 and eventually hits it, after which one obtains a pair of singularities ?7r/3 it a of which 
the right-moving one survives in the kink limit. Due to ( [E.3C| ) this leads to a zero of the 
Y{9) sX 9 = a, which means that the pseudo-energy function e{9) diverges to — oo at the 
same point. To treat this problem we chose the prescription given in and redefined the 
pseudo-energy by extracting the singular term: 

ei{e) = ti{9) - logtanh^(e' - a) 



^To verify that the solution is correct one can compare it to numerical results from TCSA ]lE 
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so that the integral equation for the massless hmit becomes 

ei(^) = - ^^{9 - e') log {tanh-^ie' - a) + e'^'^^^')) 

It is supplemented by the following condition for the singularity position 
/- „ 3 /■°° d^' cosh 2(^' - a) , / ,3,,, , . 

where f denotes the principal value of the singular integral. These equations can now be 
iterated as described above and Yi can be expressed as 

Yi{e) =tanh|(0-a)e^"i(^) 

The singularity position for this state is 

a = 0.49577315 . . . 
and the value of the effective central charge is 

ceff = ^e" + dee'\og\l + Y,{er'\ = 

TT vr^ 7_oo 5 

(due to Yi{a) = 0, the integrand has a logarithmic singularity, which is integrable, but makes 
a precise numerical evaluation harder). 

For the 4th excited state corresponding to the energy momentum tensor T = L_2id, I4 
has a real zero a together with an ordinary complex 6, so that 

3 



Yi{e) =tanh^(0-a)e^"4(^) 



where 64 satisfies 



=e' + log -^-^ - ^^{6 - 9') log (tanh ^(^' - a) + e--^(^')) 
The positions of the singularities can be determined from the equations 



TT = V3e"-i 



, S{a + m/?>-e) , S{a-iTTl?>-e) 
log 777 T— ^ - log 



S{a + iT:/?,-e*) S{a-i7T/3- 9*) 

3 d9' cosh2{e' - a) ^ ( 3,,, , 
^-.L ^ smh3(V-a) (*^"^ 4 " + ' 

= i + n(^) 

with the results 

a = 0.83331625 ... 9 = 2.49331757 • • • + 0.52367139 ■■■ xi 
and the effective central charge is 

6\/3 „ 12i / ^ , 6 



Ceff 



TT 

262 
'~5~ 



e" + 1^ fe^" - e^"*) + % T d9e' log |l + Y^-^ 
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